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Abstract. We present a set of three unified equations of states (EoSs) based on the 
nuclear energy-density functional (EDF) theory. These EoSs are based on generalized 
Skyrme forces fitted to essentially all experimental atomic mass data and constrained to 
reproduce various properties of infinite nuclear matter as obtained from many-body cal- 
culations using realistic two- and three-body interactions. The structure of cold isolated 
neutron stars is discussed in connection with some astrophysical observations. 

Of 



1. Introduction 

Neutron stars, born from the catastrophic gravitational collapse of massive stars with 
M > SM Q at the end point of their evolution, are among the most compact objects 
in the Universe. The extreme conditions encountered in their interior make their de- 
scription a very challenging task. Indeed, neutron stars are expected to contain very 
different phases of matter, from ordinary nuclei to homogeneous nuclear matter; their 
core with densities exceeding several times the nuclear matter saturation density might 



contain additional particles like hyperons or even deconfined quarks (Ha ensel et al. 



20071) . The EDF theory allows for a consistent and computationally tractable treat- 
ment of these various phases. We have determined the global structure of neutron stars 
using three different unified EoS s based on the recently developed EDFs BSkl9, BSk20 



and BSk21 (iGorielv et al.ll2010h . 



2. Unified equations of state for cold neutron stars 

The interior of a neutron star, assumed to be in full thermodynamic equilibrium at zero 
temperature, is composed of three main regions: (i) the outer crust (at densities above 
a few times 10 4 g cm -3 and below ~ 4 x 10 11 g cm -3 ) where nuclei arranged in a 
body-centered lattice coexist with a quantum gas of electrons, (ii) the inner crust where 
neutron-proton clusters coexist with both electrons and unbound neutrons and (iii) the 
liquid core (above ~ 10 14 g cm" 3 ) which consists of a uni form mixture of nuc leons 



and leptons. No other particles will be considered here (see iChamel et al.H2012l for a 
discussion of the consequence of a possible phase transition in the core). 

These three distinct regions ar e described consiste ntly using the generalized Skyrme 
EDFs BSkl9, BSk20 and BSk21 (IGorielv et all 1201 0l) . These EDFs were fitted to the 
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2149 meas ured masses of nu clei with N and Z > 8 given in the 2003 Atomic Mass 
Evaluation dAudi et al.l 12003 ). with an rms deviation of 0.58 MeV. The masses were 
obtained by adding to the Hartree-Fock-Bogoliubov (HFB) energy a phenomenological 
Wigner term and correction term for the collective energy. The EDFs were also con- 
strained to reproduce various properties of homogeneous nuclear matter as obtained 
from many-body calculations using realistic two- and three- nucleon interactions. In 
particular, the three different EDFs BSkl9, BSk20 and BSk21, were fitted to three dif- 
ferent neutron matter EoSs, reflecting the current lack of knowledge of the high-density 
behaviour of dense matter. More specifically, BSkl9 (BSk21) was adjusted to the soft- 
est (stiffest) EoS of neutron matter known to us, whereas BSk20 was fitted to an EoS 
with an intermediate stiffness. 

T he EoS of the out er crust has been determined in the framework of the BPS 
model dBaym et al.lll97ll) using the latest experimental atomic m asses complemented 
with theoretical masses obtained from our HFB mass model (see IPearson et al.ll201 1 . 
for details). For the inner crust, we have applied the Thomas-Fermi method extended 
up to the 4th order with proton she l l corrections a dded pertubatively via the Strutinsky 
integral theorem dOnsi et al.ll2008l : IPearson et al.ll2012h . This method is a high-speed 
approximation to the self-consistent Hartree-Fock equations. Neutr on shell corrections , 
which are known to be much smaller than proton shell corrections dChamel et al.l2 007). 
have been neglected. In order to further reduce the computational work, nuclear clus- 
ters have been assumed to be spherical and the Wigner-Seitz approximation has been 
adopted to calculate the Coulomb energy. At high densities, these clusters are found 
to dissolve into a plasma of neutrons, protons, electrons and muons thus delimiting the 
boundary between the inner crust and the core. 



3. Neutron star structure 

In this Section, we discuss the global structure of cold isolated neutron stars in con- 
nection with some astrophysical observations. For this purpose, we have solved the 
general-relativistic Tolman-Oppenheimer-Volkoff equations with our unified EoSs BSkl9- 
20-21 using the LORENE librarjQ. 

The gravitational redshift z sur f = (1 - 2GM/Rc 2 )~ 1 ^ 2 - 1, where M is the grav- 
itational mass and R the circumferential radius of the star, can be obtained from the 
identification of spectral features in the electromagnetic radiation from the surface of 
neutron stars. Fig. Q] shows the surface redshift as a function of M for our three EoSs 
BSkl 9-20-21. For compari son, we have also plotted the redshift predicted by the SLy4 
EoS (Chabana t et al.|[l997 ). The shaded area is prohibited by General Relativity and 



causality. The horizontal band bounded by dashed lines corresponds to the redshift es- 
timated fr om the e + e~ annihilation line in the spectrum of the gamma-ray burst GRB 
790305b dHigdon & Lingenfelterill990l) . Finally, t he horizontal dashed -dotted line at 
Zsurf = 0.35 corresponds to the redshift obtained by ICottam et all (2002) from absorp- 



tion lines in the X-ray spectra of the low mass X-ray binary EXO 0748-676. The 
redshifts predicted by our EoSs are compatible with these observations but imply dif- 
ferent neutron-star masses: the stiffer the EoS is, the larger is the mass for a given 
redshift. 



*http://www.lorene. obspm.fr 
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Figure 1 . Gravitational surface redshift as a function of the gravitational mass for 
the EoSs BSkl9-20-21 and SLy4. 



Figure|2]shows the mass-radius relation for non-rotating stars. The horizontal band 
around 2M Q indicat es the recently measu red value of the mass of PSR J1614-2230 
including error bars (|Demorest et alJl2010h . The maximum mass of neutron stars pre- 
dicted by our EoS BSkl9 is fou nd to be too low, e ven after taking into account the 
effects due to the pulsar rotation (Cham el et alJl201lh . 



4. Conclusions 



We have developed a set of three different unified EoSs of cold isolated neutron stars 
in the framework of the EDF theory. These EDFs are all based on generalized Skyrme 
forces, fitted to essenti ally all experimenta l nuclear mass data with an rms deviation 
falling below 0.6 MeV (lGoriely et alj|2010h . At the same time, these EDFs were con- 
strained to reproduce several properties of homogeneous nuclear matter (including the 
neutron-matter EoS) as obtained from microscopic calculations using realistic nucleon- 
nucleon potentials. 

These EoSs have been applied to compute the global structure of neutron stars. 
Even though our softest EoS BSkl9 seems to be favor ed by measure ments of the n~ /n + 
production ratio in heavy-ion collision experiment s (IXiao et al.ll2009i). it is ruled out by 
the recently measured mass of PSR J1614-2230 dDemorest et alJl201Ch . This conclu- 
sion mi ght suggest that the dense core of neutron stars could be made of non-nucleonic 
matter (Cham el et al .1120121) . 
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Figure 2. Gravitational mass as a functional of the circumferential radius for the 
EoSs BSkl9-20-21 and SLy4. 
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